Kondo temperature for a quantum dot in an Aharanov— Bohm ring 
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We study the Kondo temperature of a quantum dot embedded into one arm of an Aharonov- 
Bohm interferometer. The topology of a disordered or chaotic Aharanov-Bohm ring leads to a 
stochastic term in the scaling equation and in the renormalization procedure. As a result, the 
Kondo temperature displays significant fiuctuations as a function of magnetic flux. 

PACS numbers; 72.15.Qm, 73.63.Kv, 73.23.Hk 



X3 
O 

o 



> 
m 

o 
o 



I 

o 
o 



X 
S3 



Recently, the transmission phase of electrons passing 
through a quantum dot (QD) embedded in an Aharanov- 
Bohm (AB) interferometer has been measured 0. This 
was done in the Kondo regime and over an energy in- 
terval given by the width of the "odd" valley separating 
two conductance Coulomb-blockade peaks. The QD with 
half-integer spin plays the role of the Kondo impurity . 
Theory predicts the phase shift to be tt at low tempera- 
ture The results are puzzling and do not follow this 
prediction. In the present Letter we show that the T- 
matrix 7k of the Kondo problem acquires specific prop- 
erties due to the topology of the AB ring. These have not 
been taken into account previously and profoundly affect 
the dependence of 7k on the AB phase. They serve as 
one example for the role of finite-size effects and meso- 
scopic fluctuations which, together with experiments 
on semiconductor quantum dots ^\ and nanotubes 
have rekindled interest in the Kondo effect recently. 
Implications of our results for the experimental analysis 
are briefly discussed. 

Model. We consider a QD carrying an odd number of 
electrons in the Kondo regime. The QD is embedded into 
one arm of an AB ring threaded by a weak magnetic flux 
$. We take account only of the AB phase </> = 27r$/$o- 
Here $o is the elementary flux quantum. The AB ring 
is attached to A leads labelled j — 1, . . . , A. Without 
counting spin degrees of freedom labelled s = ±1/2, lead 
J carries Nj channels (transverse modes with energies 
below the Fermi surface) labelled a — l,...,Nj. The 



total number of channels is N 



The leads are separated from the interior of the AB 
ring by fictitious barriers, and the interior from the 
QD by real tunneling barriers caused by the gate po- 
tentials. With E the energy, the lead Hamiltonian is 



L = EU r ^'^.^ ^ ^ ^]asEC,asE, with 

and c the creation and annihilation operators obeying 
{Cja.E'S'a'.'iS'} = 5 j j, 5aa' 5 ss' 5{E - E') . The Hamilto- 
nian -ffAB of the AB ring has a discrete spectrum with 
eigenvalues with /i — 1,2,... and is given by i?AB 
where {d\,,d^^,} = 5. 
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'(5^1/. Leads 



and AB ring are coupled by the tunneling Hamiltonian 
^^AB-L = Esj I dE W,a-AE) HasEd,. + H.c.]. We 
assume that the real tunneling matrix elements Wja-f_i {E) 



change little over an energy scale given by the width of 
the Kondo resonance. We also assume that the levels 

lie so dense (and that their coupling to the leads is 
so large) that all scattering through the AB ring in the 
absence of the QD is smooth in energy on the same scale. 

The tunneling matrix elements vj^ between the QD 
and state /i of the AB ring carry an upper index P = 
L, R denoting tunneling through the left (right) barrier, 
respectively. The occurrence of two amplitudes for each 
state fi of the AB ring refiects the topology of the ring. 
We choose the matrix elements vj^ real and display the 
AB phase explicitly in the Kondo Hamiltonian. Gauge 
invariance allows us to put the entire effect of the AB 
phase 4> onto the right barrier separating the QD from 
the AB ring. Whenever the electron leaves (enters) the 
QD through that barrier, it picks up the phase factor e"^ 
(e"*"^, respectively). We define v^=v^ -\- e^'^v^. 

We focus attention on the Coulomb blockade mid- 
valley region. Then it is justified to use the s-d model 
for the Kondo resonance 8]. (We have convinced our- 
selves that our arguments hold likewise for the Anderson 
model 0- However, we have not yet performed a com- 
plete analysis of the latter model) . In the s-d model, a 
spin S (which represents the total spin of the electrons 
on the QD) is coupled to the spins of the electrons in the 
AB ring. 



Hk = 



Er 



■^Vf,dl^{ass' ■ S)v*d^ 



(1) 



Here Ec is the charging energy of the QD, and a stands 
for the vector of the three Pauli spin matrices. The total 
Hamiltonian is the sum of these terms, 

n = H^ + HAB + Hab-l + Hk . (2) 

Using standard techniques we obtain for the non- 
diagonal dimensionless conductance coefficients gji with 
j ^ I oi the many-body Hamiltonian Ti 



9ji 



AE 



AE j 



T,[E)G^{E)Vi{E)G''{E) . (3) 



Here / is the Fermi function, C and G" are the ordinary 
retarded and advanced equilibrium many-body Green's 
functions of TL, and [Fj]^^ = 2i:Y,a Wja:tj.Wja-^. 
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In the absence of the Kondo term, we deal with a pure 
single-particle problem, with a Green's function Ggl^^, 
given by [Gl\-^ = {E~e^)5^,^ + {i/2)Y,j\^j\p.u, and non- 
diagonal elements of the spin-independent scattering ma- 
trix S%{E) = -2iT:Y.^.W,aAGlUWib,. [11 [HI. 
Hence, as expected, the conductance is given by the 
Landauer-Biittiker formula. 

We return to the many-body case and expand G^^ [E) 
in powers of to write G^^^{E) = Gq.^^,{E) + 
[G5(ii;) rKGS(i?)] where 



TK = H^Y.^Gl{E)H^r 



(4) 



defines the T-matrix for the Kondo effect. Using 7k and 
the definition of S'^^\ we write gji as 



93l = 



AE 



dE 



Tr. 



ab I 



c.c 



-Tr, 



V Wja-AE)[Gl Tk Gl]^,WibJ 



(5) 



Three processes contribute to transport through the AB 
ring, each represented by one term on the r.h.s. of Eq. |(SJl: 
Scattering of an electron through the AB ring without 
its passing through the QD, the interference term be- 
tween the amplitude for that process and the amplitude 
for (multiple) visits of the QD, and the square of the lat- 
ter amplitude. This structure differs fundamentally from 
that for a QD in the Kondo regime connected to two 
leads without the topology of the AB ring. Here only the 
interference term occurs [3. This is caused by the spe- 
cial form of the coupling between Kondo resonance and 
channels. 

Dependence upon the AB Phase. We compare the de- 
pendence of gji on the AB phase with that obtained 
for a simpler problem: An AB ring with a QD that 
carries a single level at energy Eq and causes a Breit- 
Wigner resonance The resulting expression for g^i 

(j ^l) has the form of our Eq. (|3J| except for the replace- 
ment im -2mWja-AE) [Gl{E) Tk Gl{E)]^, Wib,, ^ 
-i{lja + rijae"l') {E-Eo+iT/2)-^ {jw + Vibe-'^) ■ The 
amplitudes jja and ijja are complex and independent of 
energy E and AB phase 0. The only dependence of the 
resulting expression on the phase </> is that given explicitly 
and that of the total width T = J^ja \"fja. + e^'^Vjal'^- Un- 
der the assumption that the electrons move chaotically 
and/or diffusively in the AB ring, the dependence of F on 
(j) is ^3 of order 1/\/N and, thus, negligible whenever 
TV ^ 1. The experimental determination of the trans- 
mission phase through the QD (defined as the phase of 
the Breit-Wigner term) is, thus, possible for ^ 1. 



In the case of Eq. jSJ , the T-matrix 7k depends non- 
trivially on 0, see Eqs. Q and (0. Thus, the central 
question is: How strong is this dependence of 7k on (jf! 
The theoretical treatment of the Kondo resonance re- 
quires scaling and/or renormalization techniques. Our 
question can only be answered in this framework. 

Poor Man's Scaling. The essential aspects of poor 
man's scaling emerge already when one considers the ma- 
trix element of the T-matrix to second order in , given 
by {vs' ITJkGJK-E) HK\fJ.s). We suppress the external fac- 
tors Vf^ and focus on the intermediate-state summation 
involving Gq{E). Two processes contribute to this term: 
The scattering of an electron and that of a hole. Aside 
from spin-dependent factors, the first process yields 



Tp T 

Vr ^_ [- I Wr;jaFja;lbWlb,p 



(6) 



E+- 



where [F-^]ja;ib ^ SjiSab [E+ - e^]-^W^,ib. 

A similar decomposition into a term involving only the 
states in the AB ring and a remainder applies to hole 
scattering. Our problem differs from the standard Kondo 
situation of a magnetic impurity coupled to the conduc- 
tion band in three main respects, the last two being 
related to the topology of the AB ring: (i) The energies 
Er do not lie in the conduction band but are defined as 
the eigenvalues of electrons moving in the AB ring. Thus, 
the summation over r extends to +00. (ii) Gq{E) differs 
from the Green's function for free electrons. This leads 
to the 2"^^ term on the r.h.s. of Eq. ©. (iii) The w^'s 
involve the AB phase (j). We address these in turn. 

(i) The coupling matrix elements and v^- reflect the 
dynamics of the coupling of the QD to the AB ring. The 
Vr^s are overlap integrals involving the wave function x 
of the electron on the QD and that of states r with mean 
level spacing A in the AB ring. Two tendencies char- 
acterize the behavior of the Ur's. (a) As the number of 
nodes of the state r in the overlap region becomes much 
bigger or smaller than that of x, the overlap is strongly 
reduced, (b) In a diffusive and/or chaotic AB ring, the 
states r within an energy interval given by the Thou- 
less energy Et, are strongly mixed. We take account of 
both features in terms of a cutoff model. Within an en- 
ergy window of bandwidth D < Et around the energy 
of the state Xj the Wt's are uncorrelated Gaussian 
distributed random variables with mean value zero and 
a common constant variance while the Vt's for states 
T outside of this window vanish. For both ballistic and 
diffusive rings Et is large compared to the mean level 
spacing A of the states in the AB ring but small in com- 
parison to the Fermi energy. The cutoff then applies 
likewise to the summation over particles and to that over 
holes. The cutoff restricts the summation over r to a 
finite number g « Et j ^ states, with g^\. To identify 
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the influence of the large contributions to Eq. © and to 
resum them, we let g —^ oo and use scaling. 

(ii) Scaling needs to be applied only to the first term 
on the right-hand side of Eq. © since only this term pro- 
duces a logarithmic singularity at the upper band edge. 
As for the second term, the expression in brackets is ob- 
viously finite. There are two contributions to the sum- 
mations over T and p. For t = p, we sum squares of 
amplitudes Vr but also of the propagator. The result is 
inversely proportional to D and vanishes for D oo. For 
T p, we use that the f^'s are Gaussian random variables 
with mean value zero. Their random signs quench the 
intermediate-state summations and ensure that the re- 
sult is convergent even as g ^ oo. Thus, the singular con- 
tributions to the second-order particle matrix element 
are obtained by replacing Gq{E) by [E~^ — Hab]~^- Cor- 
responding statements apply to the terms of higher order 
in Hk- These arguments apply likewise to hole scatter- 
ing. Thus we have reduced our problem to one which is 
quite similar in appearance to the standard Kondo case: 
Singular terms arise only from the coupling of the QD 
to the states in the AB ring, i.e., from the Hamiltonian 
Hab + Hk- The terms and Hab-l in H of Eq. © 
are irrelevant for scaling. 

(iii) We turn to the AB phase appearing in the 
intermediate-state summation. The first term on the 
right-hand side of Eq. © contains the phase through 

= (ui)2 +{v^)'^ -h2cos(/) v!^v^. To define the scal- 
ing variable, we write this expression in the form [{v^)'^ 
-|-(t;.^)^][l -f cosc^^cos^/^t] where 



cos ipr 



(7) 



Since the u^-'s are uncorrelated Gaussian variables, '0r is 
also a random variable which is uniformly distributed in 
the interval < f/'r < ti"- Moreover, the ipr's with differ- 
ent indices p are uncorrelated. We replace the positive 
quantity [(w^)^ -I- (w^)^] by its average 2v'^. An equiva- 
lent assumption (independence of the coupling strength 
on the index r characterizing states in the conduction 
band) is usually made in the standard Kondo case. We 
checked that this assumption does not affect the generic 
features of our results. 

We define the scaling variable J = Av^ /En and apply 
the standard algebra of poor man's scaling [Sj to write 



5J_ 
SD 



2poJ^ 
D 



(1 + cos cos V') 



(8) 



Here po — A^^ is the density of states \p) in the AB 
interferometer. We assume that po is constant. Because 
of the appearance of the random variable ip, Eq. (jS]) is 
a stochastic differential equation. To solve it, we first 
consider the case without stochastic term and define the 
"standard scaling trajectory" by choosing a Kondo tem- 
perature Tk and initial values pqDq (where Dq « Ex) for 



the bandwidth and pqJq for the coupling strength such 
that the relation Tk = Dq exp(— l/2po>/o) holds. We 
then solve Eq. ((HJ by reducing poD in unit steps drawing 
-0 from a random-number generator. We thereby elim- 
inate the states of the AB ring one by one. In the p^^ 
step, the coupling constant poJ is "kicked" out of one 
scaling trajectory and into another by the term propor- 
tional to cos^cosi/'p- The effect becomes larger as p^D 
is decreased and 6J/6D increases. As a result, every re- 
alization {ipp} of the sequence of i/j-values determines a 
different scaling trajectory given by a discrete sequence 
of points {Dp, Jp}. The scaling process has necessarily to 
end when all ring states are reduced to a single remain- 
ing one that absorbs all interaction effects, giving a lower 
energy cutoff A. The lower cutoff is modified for two 
reasons: (i) The AB ring is open, and the states p have 
finite widths with mean value F — AiV/(27r). We have 
neglected F by omitting the second term on the r.h.s. of 
Eq. jnj. This is justified only as long as » F. (ii) 
The electrons which screen the QD spin, form a Kondo 
cloud of size ^ in the AB ring. The energy scale E^ which 
is required to attain a state of correlation length ^, can 
be estimated using the uncertainty relation. Hence, the 
lower cutoff is A = max{A, F, i?^}. 

We have solved the discretized version of Eq. ((HJ for a 
choice of parameters which resembles the weak Kondo 
coupling experiment in Ref. Q: (i) F/A w 10 (A = 
6,iVj « 10); E^/A « 1 (since ^ w 1/im is about the 
size L of the ring), (n) Et/A = ^/AttT^ where is the 
number of electrons on the ring. This formula follows for 
ballistic system from Et h/L and the Weyl formula. 
The same formula yields « 1000. (in) Tk/A « 5. 

Our results are shown in Figure ^ The insert shows 
five representative trajectories for cos (f> — 1 (where 
the stochastic effect is largest). The deviations from 
the standard scaling trajectory increase with decreasing 
bandwidth, as expected. Averaging over the stochastic 
trajectories, we recover the standard scaling trajectory 
(straight line in the insert). In the main part of the 
figure, we show the distributions P{Tk) of Kondo tem- 
peratures Tk for cos(/) =1/2 and coscf} = 1. These are 
obtained from 10"* realizations of {V'p}- For each realiza- 
tion, the final value of Tk is found by a best-fit proce- 
dure. Mean value {Tk) and dispersion 5Tk are deduced 
from P{Tk)- The average value {Tk) is very close to the 
standard scaling value and does not show a significant 
dependence on cos </>. The fluctuations are large, how- 
ever, and 5Tk increases with cos 0, reaching a maximum 
value of STk/{Tk) « 0.25 at cos(/> = 1. 

The numerical simulations show that the fluctuations 
become larger (smaller) as either pqTk or the lower cut- 
off A are decreased (increased). The reason is that the 
stochastic "kicks" become more effective for smaller val- 
ues of Tk and/or D, see Eq. Our results are less sen- 
sitive to Et- As Et is increased, the slopes of the scaling 
trajectories become dominated by the upper band edge. 
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FIG. 1: Distribution of Kondo temperatures Tk, for cos(f) = 
0.5 (narrow distribution) and cos0 = 1 (wide distribution), 
both for (Tk)/A — 5. Each histogram contains 10* scahng 
trajectories. Insert: Typical scaling trajectories for cos = 1. 

were the stochastic kicks have little effect. However, this 
happens very slowly: We observe only small quantitative 
changes as we increase Et by an order of magnitude. 

For the physical interpretation of our result, we say 
that an electron is on the QD when it interacts with the 
fixed spin S of the Kondo Hamiltonian iJx- Each term 
of the perturbation expansion describes a process or a 
sequence of processes where an electron leaves the QD, 
moves in the AB ring, and re-enters the QD. Whenever 
the electron leaves and enters the QD through the same 
barrier, the result is a factor (v^Y oi' {'^^)'^ i as the case 
may be. Processes of this type are not different from 
the ones that are taken into account in the standard ap- 
proach and lead to the standard form of the scaling equa- 
tion. The fact that the QD is coupled to the AB ring via 
two barriers only contributes a factor two and is other- 
wise irrelevant. There are, however, additional processes 
where the electron leaves the QD through the left bar- 
rier and re-enters it through the right one, or vice versa. 
Such processes correspond to a complete loop through 
the AB ring, carry the factor v^v^ and the AB phase, 
and are absent in the standard approach. Loops in either 
direction are equally likely, hence the dependence of the 
relevant term on cos0 v^v^. The AB ring is assumed to 
be chaotic or diffusive. This makes the loop terms in the 
scaling equation stochastic. 

Renormalization. Poor man's scaling breaks down for 
temperatures below the Kondo temperature Tk- We dis- 
cuss this temperature regime only qualitatively in order 
to demonstrate the existence and action of the stochastic 
term. We use renormalization group arguments with- 
out actually performing the calculation. Following stan- 
dard procedure one would construct a new basis 
of orthonormal states |n) where n = 0,1,2,.... With 



|vac) denoting the vacuum, the state |0) is defined as 
|0) = (E.k.^n"'EM<4|vac). For n > 1, the states 
\n) would be obtained with the help of Schmidt's or- 
thogonalization procedure from the sequence Hxb\0), 

i/|^g|0), (We have shown that it is again justified 

to approximate the Hamiltonian Ti. — Hk by Hab)- The 
terms in the resulting tridiagonal Hamiltonian |3| depend 
upon the matrix elements via expressions of the type 
J2f^ w^(e^)™, with TO = 0, 1, 2, . . .. Just as for poor man's 
scaling, the averages (fluctuations) of such expressions 
are independent of (dependent on) the AB phase. The 
fluctuations are maximal for cos (p — 1 and vanish for 
cos (f> = 0. Again, the stochastic part is expected to cause 
a spread in the renormalization group trajectories. The 
trend of the spread with temperature T can be ascer- 
tained as follows. The Kondo temperature Tk defines a 
marginal fixed point of the problem, while T = is a sta- 
ble fixed point. All trajectories converge to that point. 
Therefore, as the sample is cooled down, the stochastic 
spread of the trajectories first increases monotonically to 
reach its maximum at T = Tk. Thereafter, the spread 
decreases monotonically and vanishes at T = 0. 

Discussion. Because of the topology of the AB ring, 
the Kondo temperature Tk of a QD embedded into the 
ring depends stochastically on the AB phase (j>. The re- 
sulting spread in Kondo temperature versus AB phase 
is significant, amounting to 20 or even 40 percent. The 
influence of stochasticity is expected to decrease as the 
temperature is lowered below Tk, and to vanish at T = 0. 
The precise form of the dependence of Tk on can- 
not be predicted for an individual experimental sample. 
We expect that for T « Tk, this stochastic dependence 
seriously limits attempts to determine the transmission 
phase through the QD experimentally, although quanti- 
tative predictions will have to be based upon a study of 
the Anderson model. Such work is now in progress. 

This work was partly done while one of us (HAW) vis- 
ited Brazil under the auspices of a CAPES - Humboldt 
award. CHL is supported by CNPq and Instituto do 
Milenio de Nanociencias. HAW acknowledges stimulat- 
ing discussions with Y. Gefen, Y. Imry, and Y. Oreg. 
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